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Abstract: The possibility of diquark condensation at sufficiently large baryon chemical 
potential and zero temperature is analyzed in QCD at strong coupling. In agreement with 
other strong coupling analysis, it is found that a first order phase transition separates a 
low density phase with chiral symmetry spontaneously broken from a high density phase 
where chiral symmetry is restored. In none of the phases diquark condensation takes place 
as an equilibrium state, but, for any value of the chemical potential, there is a metastable 
state characterized by a non-vanishing diquark condensate. The energy difference between 
this metastable state and the equilibrium state decreases with the chemical potential and 
is minimum in the high density phase. The results indicate that there is attraction in the 
quark-quark sector also at strong coupling, and that the attraction is more effective at high 
baryon density, but for infinite coupling it is not enough to produce diquark condensation. 
It is argued that the absence of diquark condensation is not a peculiarity of the strong 
coupling limit, but persists at sufficiently large finite couplings. 
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1. Introduction 

Hadron phenomenology suggests that quarks confined within the baryons are strongly 
clustered into diquarks ||, [|]. Theoretically, a weak coupling analysis of QCD at the 
level of one gluon exchange shows that the quark-quark scattering amplitude is attractive 
in the color anti-triplet channel and repulsive in the symmetric color sextet channel. The 
attraction in the color anti-triplet channel due to one gluon exchange provides a mechanism 
for bounding the diquark in a state the color of which is neutralized by the remaining quark 
within the baryon. The quark-quark interaction induced by non-perturbative instanton 
effects produces also attraction in the color anti-triplet channel and repulsion in the color 
sextet channel. 

The effect of the quark-quark attraction should be more dramatic at high baryon 
densities ||. At zero density the loosely bound diquarks are in its turn bound to quarks 
producing the color singlet baryons. At high baryon density, however, the color charge is 
expected to be deconfined and the weakly interacting quarks may form a Fermi surface 
that the quark-quark attraction will render unstable. Then, the phenomenon of color 
superconductivity would take place via a BCS mechanism, with the formation of a diquark 
condensate f|, ||]. The diquark condensate can be studied in QCD in the weak coupling 
regime via the Schwinger-Dyson equations pL ^, ||, |9| and instabilities of vertex functions 
[|To|] . The diquark condensation offers the possibility of exotic phenomena at high baryon 
density, as color-flavor locking and unlocking and crystalline superconductivity |11, 12, 13 



14], and the continuity of the quark and hadronic phases Jl5| ]. 

The diquark operator carries color charge and is not gauge invariant. Hence, general 



arguments forbide diquark condensation in the naive way jL(|] (see also [17]). Diquark 
condensation must be understood not as the spontaneous breaking of the gauge symmetry, 
but as a kind of Higgs mechanism that has the observable consequencies of color supercon- 
ductivity ||,|||. 
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Recently, strong evidence of diquark condensation in two colors QCD has been found 



|2C, 21, 22, 23, 24, 25]. The analysis in this case is simpler than in three colors QCD for 
two reasons: there is a gauge invariant diquark operator and the fermionic determinant is 
positive at finite chemical potential, so that Monte Calo simulations are feasible. 

In three colors QCD, however, the arguments in favour of diquark condensation at high 
density are based on weak coupling analysis (although nonperturbative effects are taken 
into account via instantons), that should be valid at asymptotically large densities. It is 
very important to verify the diquark condensation relaxing the weak coupling assumptions. 
Unfortunately, numerical simulations of lattice QCD at finite density are unfeasible ought 
to the sign problem. Hence, the strong coupling techniques are valuable to get insight in 
the problem. Recent attempts to study QCD at finite density in the strong coupling limit 



by using Hamiltonian techniques have been developed in [26, £7], 28]. In this paper, we 
shall study the possibility of diquark condensation in the strong coupling limit using the 
path integral formalism. 

2. Effective action at strong coupling 

Let us consider QCD regularized on a euclidean four dimensional lattice with one flavor 
of staggered fermions. The gauge group is SU(3). Let us define the meson and baryon 
operators: 

M{x) = i> a (x)i> a (x) (2.1) 

B(x) = \e abc r{x)^\x)^ c {x) (2.2) 
b 

B{x) = \e ab 3 a (x)^ b (x)r{x) (2.3) 
6 

where summation over repeated color indexes, a, b, c, is understood. At strong coupling, 
1/g 2 = 0, the gauge field can be integrated out and gives the following effective action for 



the fermions [29]: 



S s f = m J2 M ( x ) ~ \ E M (^m(x, y)M{y) - £ B(x)V B (x, y)B{y) 

x x,y x,y 

+ 5^6 M 2 (x)M 2 (x + u) - B{x)B{x)B(x + u)B{x + v) , (2.4) 

x,v x,v 

where mo is the bare fermion mass and 

Vu(x,y) = — y^X S v,x+" + ( 2 - 5 ) 
Vb{x,v) = j^I/i/WW _ fA x )~ ls y,x-u] , (2.6) 

V 

where v are the unit vectors in the four space time directions and 

f e 3 ^ v = n 
r)v{x), v = 1,2,3 
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where [i is the baryon chemical potential and r\ v {x) are the Kogut-Susskind phases. 
The zero temperature grand canonical partition function is given by 

Z = f [dijjdiP] exp[-S sc ] ■ (2.8) 

We will neglect the last two terms of the fermion effective action, which involve eight- 
fermion and twelve-fermion vertices. It has been proved that these terms are sub-dominant 
in a 1/d expansion, where d is the space-time dimension [^9| . In d = 4 their contribution to 
the thermodynamics is small and does not change qualitatively the phase diagram |3l| . We 
linearize the remaining meson and baryon terms in the usual way, introducing a bosonic 
field, (f(x), and two fermion fields, b x and b x , via the Hubbard-Stratonovich transforma- 
tions: 

exp{^(M,V M M)} = (dety M )" 1/2 y[^]exp{-i(^,y M V) - (<P,M)}, (2.9) 

exp{(5, V B B)} = det V B J [dbdb] exp{-^(6, V^b) + (b, B) + (B, b)} . (2.10) 

where we use the notation (f,Ag) = 52 xy f(x)A(x,y)g(y) and (f,g) = J2 X f(x)g(x). 
Let us introduce the diquark fields 

D a (x) = e abc ^ b {xW{x) (2.11) 
D a {x) = e ab 3 b {xW{x). (2.12) 

Note that D a {x) transforms as a quark (color triplet) and D a (x) as an antiquark (color 
anti-triplet) under gauge transformations. It is convenient to our purpose to rewrite the 
term exp[(6, B) + (B, b)] entering equation ( p.lOj) by using the following identity: 

exp[b x B(x) + B(x)b x ] = exp[^b x b x M(x) + 2M 2 (x)] J # f (x)#(x)exp \-<j>\(x)(i> a {x) 

-<t>l(x)[~b x r(x) + D a (x)] - M^)[\r{x)b x + D a (x)]\ ,(2.13) 

where cp a (x) is a complex bosonic field that carries color charge in the fundamental rep- 
resentation. The above series of transformations lead to the following expression for the 
grand canonical partition function: 

Z = J [d$d$][&db][d*p][d<ltd4>] exp[-S{jf] , (2.14) 

where 

Sfb =^2<PUx)Mx) + \^2^(x)V M 1 (x,yMy) + ^2b x V B \x,y)b y - lndetV B 

x x,y x,y 



+ Yl { ( m + ^ " ^Kb x )M(x) - 2M 2 (x) 

■T 



+ 4 > \{ x )[h x r{x)+D a {x)\+U^)[\r{x)b x + D a {x)} ) . (2.15) 
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Notice that the scalar field 4> a (x) acts as a Higgs field with color charge in the fundamental 
representation. 

The integral over the fermion fields ip a (x) and ip a {x) factorizes as a product of integrals 
at each lattice site and can be readily performed. Afterwards, the auxiliary fields b x and b x 
can be integrated out, and we arrive to a representation of the grand canonical partition 
function as a functional integral over bosonic fields: 

Z = J [d<p}[d(f)U<t)] exp[-Sf] , (2.16) 

with 

Sf = 5>l(*)^a(s) + \j2<p(x)Vu\x,y)<p(v) - bidets - V^Q 2 } , (2.17) 

x x,y 

where &i(x,y) = ®i(x)5 xy , i = 1,2, with 

6x(x) = m(x)[m 2 (x) - 4|^(x)| 2 + 12] , (2.18) 
®2(x) = ^-[4\0(x)\ 4 + m 2 (x)\0(x)\ 2 + 8|0(x)| 2 - 3m 2 (x) - 12] . (2.19) 

In the above equation m(x) = uiq + (p(x). 

The field (p(x) is related to the chiral condensate in the usual way, by 

(M(x)) = c(<p(x)) , (2.20) 

where c = (1/V) ^2 xy l^ 1 (x, y), with V the lattice volume. Hence, a non vanishing vacuum 
expectation value of <p(x) implies spontaneous chiral symmetry breaking. 

It is not completely obvious how the field 4> a {x) is related to the diquark operator. 
One way to see it is to couple sources to the diquark field and relate the derivatives of the 
partition function with respect to the diquark sources to (f> a (x). 

3. Sources for the diquark field 



The diquark fields defined by equations ( 2.1l| ) and ( [2.12 ) are not gauge invariant. Adding 



sources for these fields will make sense once a gauge is fixed. To get the effective action in 
the presence of sources, however, we do not need to fix the gauge. It can be fixed afterwords, 
when using the effective action to make computations. Let us add to the action ( |2.4[) a 
source term of the form 

J2iUx)D a (x) + 4(x)D a (x)]. (3.1) 

X 

With the same algebraic manipulations of the previous section we find that the only effect 
of the sources is to replace the functions Q\(x) and 02 (x) entering equation ( |2.17 ) by 



ej J) (x) = 8i(z) + Am(x)[4>i{x)J a {x) + Jt{x)4>a{x) - 4(x)J a (x)} , (3.2) 
G^ J) (x) = @ 2 {x) - I U\{x)J a {x) + 4(x)Ux) ~ 4{x)J a (x) 



9 



+ 



(x)\ 2 [<pl(x)J a (x) + 4(x)cf> a (x)} - <j>{{x)J a {x) 4(x)4> b (x)} . (3.3) 
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Taking the derivative respect to J a (x) and letting the sources vanish, we get the diquark 
operator expressed in terms of the bosonic fields ip and (f> a : 

D a (x) = [Mx)A(x) - i(l - \Hx)\ 2 )B(x)] 4>l(x) , (3.4) 

where 

In the above expressions matrix multiplication and inversion is understood. 
We see, as expected, that the diquark field is proportional to 4>\. 



4. Phase diagram 



The effective potential, U e g, to the lowest order semi-classical expansion (mean field ap- 
proximation) is given by the effective action ( [2.17 ) evaluated at constant fields, ip(x) = (p 



''. In this case the determinant entering equation ( 2.17j ) can be readily 



and |0(x)f = 

computed by using the Fourier transformation and we have 

1 f d 4 k 



cff 



V 2 + \p 2 



Hei-QiR'(k)], 



where 



and 



R 2 (k) 



1 

32 



e 2 



3 

£ 

z/=i 



sin 2 (^/2) -2 + e 6/ V feo +e 



-6/n — ifco 



rh(m 2 - 4v 2 + 12) 
^-[4v 4 + v 2 fh 2 + 8v 2 



3m 2 - 12] , 



(4.1) 
(4.2) 

(4.3) 
(4.4) 



with fh = too + A non vanishing expectation value, v, of the Higgs field will signal, as 
usual, the presence of the Higgs phenomenon, which, in this context, implies the diquark 
condensation. 

For the remaining of the paper we will restrict our analysis to the chiral limit, mo = 0. 
Let us analyze the phase structure as a function of the chemical potential. The effective 
potential (4.1) is invariant under the change (p — > —(p (and notice that, by definition, v > 0). 



The effective potential (4.1) has the following local minima 

,V2 



v 

V 
V 



0, 
0, 



± (V33-5) 





V ji 

V fi 

for n > 0.4416 



(4.5) 



which are independent of \x. 



- 5 - 



There is a first order phase transition at fj, c ~ 1.557. For fi < \i c the (degenerate ought 
to chiral symmetry) absolute minima are at v = and (p = ±(\/33 — 5) 1 / 2 , corresponding 
to a phase with chiral symmetry spontaneously broken. For fi > [i c the absolute minimum 
is at v = and <p = 0; chiral symmetry is restored in this phase. Obviously, the transition 
is first order. For fi < fi c the baryon density is zero, and three for fj, > fi c . Three quarks per 
lattice site is the maximum allowed by Pauli's principle, so that the transition separates 
a phase of zero baryon density from a phase saturated of quarks. The system at any 
intermediate density is thermodynamically unstable and splits into domains of zero and 
saturated baryon density. This behavior is an artifact of the strong coupling limit and has 
been observed in other strong coupling analysis at zero temperature and finite chemical 



potential (28, 32|, B|, 34|. 



The minimum at v = v5 and (p = is metastable in both phases. It describes a 
chiral symmetric state with a diquark condensate. The baryon density that correspond 
to this metastable state is a smooth function of the chiral condensate (see figure |]) that 
interpolates between zero density and saturation. This means that at any baryon density 
a state with diquark condensation can be formed and have some short life until the system 
splits into its zero density and saturated domains. The presence of the metastable state 
signals the attraction in the quark-quark channel. At strong coupling such attraction is not 
strong enough to form a stable diquark condensate. The energy difference between this 
metastable state and the equilibrium state decreases with /x, indicating that the quark- 
quark attraction becomes the more effective the higher the baryon density (see figure |2|). 
However, the energy difference is nonzero even at high \i (in the saturated phase). 

Had we computed the minimum of the effective potential with diquark sources, we 
would have obtained as metastable minimum in the limit of vanishing sources 

<t>a = Vj± . (4.6) 

Notice that for vanishing sources J a /\J\ merely selects a direction in color space. This 
does not mean that gauge symmetry is spontaneously broken, since the fluctuations of the 
redundant gauge variables will destroy the behavior of the classical level. However, if we fix 
a gauge before including the fluctuations, then the diquark field will have a non vanishing 
expectation value and the Higgs mechanism will take place. 



5. Discussion 



It is rather natural that diquarks do not condense at strong coupling, since the confining 
forces are so strong that produce almost pointlike baryons and mesons. The transition 
driven by the chemical potential separates a phase with zero baryon density from a satu- 
rated phase, with three quarks per lattice site. Hence, the system at any finite density is 
unstable and splits into domains of zero and saturated densities. Deconfinement does not 
take place at any density and diquark condensation cannot take place in a confined phase. 

It is however interesting that a metastable state in which diquark condensation takes 
place appears at any density. This shows that there is a strong attraction in the quark-quark 
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0.0 0.1 0.2 0.3 0.4 0.5 

fj. 

Figure 1: Baryon density as a function of the chemical potential in the metastable state with 
diquark condensation. 
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Figure 2: Energy difference between the diquark condensate metastable and the equilibrium state. 

channel at strong coupling too, where the weak coupling arguments based on one gluon 
exchange or instanton induced interactions do not apply. The energy difference between 
the metastable and equilibrium states decreases with the density, and is minimum in the 
saturated regime. The attraction, hence, is the more effective the highest the density. We 
expect the metastable state to have the usual features induced by diquark condensation: 
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Meissner effect, a gap in the excitation spectrum, and color superconductivity. The crucial 
question that should be addressed is whether this metastable state becomes stable at 
sufficiently large density at finite coupling. 

A partial answer is the following: assuming that the strong coupling expansion gives 
meaningful results, diquark condensation cannot take place at sufficiently large couplings, 
whatever the chemical potential. The reason is that, as we have seen, the energy difference 
between the state with diquark condensation and the stable state remains positive for any 
value of the chemical potential in the strong coupling limit. A correction to the effective 
potential of order 1/g cannot remove such energy difference if g is sufficiently large. Hence, 
diquark condensation cannot take place in the strong coupling region. If, on the other 
hand, diquarks condense at some finite chemical potential in the weak coupling regime, 
there must be a critical value of the coupling, g c , such that diquark condensation takes 
place in some interval of the baryon density for g < g c , but diquarks do not condense at 
any density for g > g c . This means that the interval of baryon densities (in lattice units) 
at which diquark condensation occurs as a stable state shrinks as g increases, and vanishes 
at g c . For g > g c the state with diquark condensate survives as a metastable state. 

It is interesting to note the different behaviour of the chiral transition: since in the 
strong coupling limit a first order phase transition separates a state with vanishing density 
and chiral symmetry spontaneously broken from a state saturated of quarks, where chiral 
symmetry is unbroken, the chiral transition at high density may persist for any value of 
the coupling, from the weak to the strong coupling regimes. The peculiarity of the strong 
coupling limit is that the chiral transition occurs at the onset chemical potential (the value 
of fx at which the baryon density starts to be non-zero). 

For two colors QCD a stable phase characterized by diquark condensation has been 
found |2^, [21], |22|, ^3|, ^4], ^5|. Two colors QCD has two peculiarities that make its inves- 
tigation simpler: 1) the diquark operator is a color singlet and can be used as an order 
parameter in the usual sense, and 2) the fermion determinant is real and there is no sign 
problem. The sign problem prevents the applicability of numerical simulations to analyze 
the phase diagram of three colors QCD and, therefore, strong coupling techniques are useful 
to improve our understanding of real QCD at finite baryon density. 
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